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Abstract —This paper proposes an encoding/decoding frame¬ 
work for achieving the optimal channel capacities of the two-user 
broadcast channel where each user (receiver) has the message 
targeted for the other user (receiver) as side information. Since 
the link qualities of the channels from the base station to the 
two users are different, their respective single-user non-broadcast 
channel capacities are also different. A goal is to simultane¬ 
ously achieve/approach the single-user non-broadcast channel 
capacities of the two users with a single broadcast transmission 
by applying network coding. This is referred to as the rate- 
diverse wireless network coding problem. For this problem, this 
paper presents a capacity-achieving framework based on linear- 
structured nested lattice codes. The significance of the proposed 
framework, besides its theoretical optimality, is that it suggests a 
general design principle for linear rate-diverse wireless network 
coding going beyond the use of lattice codes. We refer to this 
design principle as the principle of virtual single-user channels. 
Guided by this design principle, we propose two implementations 
of our encoding/decoding framework using practical linear codes 
amenable to decoding with affordable complexities: the first 
implementation is based on Low Density Lattice Codes (LDLC) 
and the second implementation is based on Bit-interleaved Coded 
Modulation (BICM). These two implementations demonstrate the 
validity and performance advantage of our framework. 

Index Terms —broadcast channel with side Information, net¬ 
work coding, nested lattice codes, LDLC, BICM 

I. Introduction 

HIS paper investigates wireless broadcast networks with 
side information at users, where a base station wants to 
deliver two different messages to two users A and B, and each 
user already has the message targeted for the other user as side 
information. Two examples in which this scenario can arise are 
as follows: 

Scenario 1: Two-way Relay Networks — In this scenario 
the message for B originates from A; and the message 
for A originates from B. These messages have been transmitted 
by users A and B to the base station via a prior uplink phase. 
The base station serves as a relay to deliver the information to 
the respective destinations, B and A, in the downlink phase. 
The current paper is related to the downlink phase in this 
setting. 

Scenario 2: Retransmission in Automatic Repeat reQuest 
(ARQ) — In this scenario, the information originates else¬ 
where (e.g., the Internet). Originally, neither A nor B had 
information on the messages for the other nodes. However, 
in prior downlink transmissions, A overheard the message for 
B, and B overheard the message for A, but neither A or B 
succeeded in decoding their own desired messages. Retrans¬ 
missions will be needed. In a future downlink retransmission, 
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both A and B have side information on the messages for the 
other user. 

Besides the above two scenarios, there are many other 
scenarios in which downlink broadcast with side informa¬ 
tion arise Network coding 0 . 0 can naturally be 
used to improve the transmission efficiency for such wireless 
broadcast networks. With network coding, the base station 
transforms the two messages into one network-coded message 
and sends the network-coded message to the two users. Each 
user then decodes its desired message by subtracting the side 
information from the received network-coded message. 

In this paper, we assume the channels from the base station 
to the two users are power-constrained additive white Gaussian 
noise (AWGN) channels. The capacities of the two channels 
can be different due to the different channel qualities from 
the base station to two users. We refer to the coding problem 
for such channels as the rate-diverse wireless network coding 
problem, because the achievable rates for the two channels 
are different; yet our goal is to be able achieve the respective 
rates with a single transmission of the network-coded message. 
Specifically, of concern are the following two questions: 

i) What is the capacity region of rate-diverse wireless 
network coding? 

ii) How to achieve (or approach) the optimal operating point 
of the capacity region using practical and implementable 
coding schemes? 

There has been prior works addressing question i). The 
capacity region of the channel under investigation has been 
identified in @-l8| using the argument of random coding; it 
was proved that the optimal point of the capacity region is the 
pair of the two point-to-point single-user channel capacities. 
However, based on the theoretical arguments in 0-0 alone, 
the answer to question ii) is not obvious. This paper is oriented 
toward question ii). We find that simple coding schemes 
can allow us to achieve/approach the point-to-point channel 
capacities, using a principle put forth by us, referred to as the 
principle of virtual single-user channels. 

To achieve the optimal point of the capacity region for rate- 
diverse wireless network coding in power-constrained AWGN 
channels, we propose an encoding/decoding framework based 
on linear nested lattice codes ||^. The merit of our framework 
is twofold. First, it shows that we can achieve the optimal 
point of the capacity region using linear structured codes. 
Second, it yields a general design principle for rate-diverse 
wireless network coding based on linear codes. We refer to this 
design principle as the principle of virtual single-user channels. 
This principle indicates that there is no need to perform joint 
two-user encoding to achieve capacities. Separate channel 
coding before network coding at the transmitter, and single- 
user decoding after network decoding at the receiver, are 
sufficient to achieve capacities. 
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Importantly, this principle is not just limited to the use of 
nested lattice codes (although the insight originally came from 
the use of nested lattice codes). For example, we can imple¬ 
ment our encoding/decoding framework using practical linear 
codes amenable to decoding with affordable complexities. In 
other words, although the framework originates from the use 
of theoretical lattice codes with inhnite dimensions, the design 
principle can be easily applied in practice using other linear 
codes. 

To illustrate, this paper shows how the optimal encod¬ 
ing/decoding framework can be applied using low density 
lattice codes (LDLC) Hg. Although the nested lattice code is 
optimal in the sense of capacity-achieving, its lattice quantiza¬ 
tion decoding is not computationally feasible. LDLC is a type 
of lattice codes amenable to low-complexity belief propagation 
(BP) decoding |TT|, Using LDLC as components, our en¬ 
coding/decoding framework for rate-diverse wireless network 
coding can be easily realized in practical systems. 

This paper further illustrates the application of the frame¬ 
work using bit-interleaved coded modulation (BICM) |T3), 
which has even lower implementation complexity than LDLC. 
BICM is already used in many practical point-to-point com¬ 
munication systems to approach capacity in the high signal- 
to-noise ratio (SNR) regime. BICM combines simple binary 
channel codes with high-order modulations to obtain different 
data rates in accordance with the SNR. This paper shows that 
applying the principle of virtual single-user channels allows 
both users to achieve their respective single-user point-to-point 
decoding performance in BICM-based rate-diverse wireless 
network coding systems. This result is contrast to prior works 
01,03’ where more complicated joint modulation schemes 
are adopted. These joint modulation schemes are designed 
for specihc modulations and are difficult to generalize. More 
importantly, as will be seen in this paper, they do not allow op¬ 
timal single-user decoding performance to be simultaneously 
achieved for both users. 

In summary, this paper has the following three contributions: 

• We put forth a framework based on structured nested 
lattice codes to achieve the optimal capacities in rate- 
diverse wireless network coding systems. This framework 
provides insight leading to a general design principle of 
using linear codes to achieve optimal rate-diverse network 
coding performance, referred to as the principle of virtual 
single-user channels. 

• We put the design principle into practice assuming the 
use of LDLC that are amenable to practical encod¬ 
ing/decoding using BP algorithms. We demonstrate the 
good performance of the resulting system. 

• We further apply the design principle assuming the use 
of BICM. The proposed BICM scheme can achieve 
the single-user performance bounds for the two users 
simultaneously. 

A. Related Works 

Ref. fTh) considered the use of linear codes to achieve the 
optimal point of the capacity region in hnite-alphabet channels 
(channels with discrete outputs). By contrast, our paper here 


focuses on power-constrained additive white Gaussian noise 
(AWGN) channels, which corresponds more closely to the 
actual physical channel at the lower layer. The coding scheme 
in G3 is a solution for rate-diverse wireless network coding 
in power-constrained AWGN channels; since it only employs 
binary linear codes and BPSK modulation, it cannot achieve 
capacity in all SNR regimes. 

The constellations proposed in p8j for noisy index coding 
can be also applied to rate-diverse wireless network coding. 
However, it aims to achieve the so-called side information gain 
but not to achieve the optimal capacity point. 

To boost the performance of rate-diverse wireless network 
coding systems, | [T4| and p?) proposed two-user joint modula¬ 
tion schemes for non-channel-coded systems. However, these 
non-systematic designs for joint modulation schemes cannot 
achieve the optimal single-user performance for the two users 
simultaneously. 

A recent independent work lig, also employed the 
idea of nested lattice codes to achieve the capacity of the 
Gaussian broadcast channel. In |Tg, 1^, a more general 
model with more than two users was considered, and a 
theoretical optimal solution was given for the general model. 
However, for the case of two users, our work and ng, 
p0| differ in the following ways. First, our encoding scheme 
employs two different coding lattices to independently encode 
the two messages, where each coding lattice is exactly the 
same as that used in traditional single-user point-to-point 
communication; the encoding scheme of 1^, on the 

other hand, jointly encodes the messages into a single coding 
lattice. Second, our decoding scheme, after the network-coding 
decoding operation, is the lattice quantization operation with 
respect to the original single-user coding lattice; the decoding 
scheme of lig, on the other hand, is a lattice quantization 
with respect to a new coding lattice representing a sub-code 
of the original lattice code. Our scheme suggests a simple 
design principle for the two-user case: the principle of virtual 
single-user channels. This principle dictates that we should 
aim for a design in which the whole decoding problem can be 
reduced to two uncoupled single-user decoding problems after 
the side information at the receivers are taken into account. 
Importantly, for a design adhering to this principle, we can 
simply use single-user channel codes, including non-lattice 
codes (e.g., LDPC codes) that are capacity approaching yet 
amenable to simple decoding. 

Overall, despite the past related works, how to achieve 
the two point-to-point single-user channel capacities simul¬ 
taneously using structured codes is not obvious for rate- 
diverse wireless network coding in power-constrained AWGN 
channels. It is the intention of this paper to hll this gap. 

The rest of this paper is organized as follows. Section II 
describes our system model. Using the argument of nestled lat¬ 
tice codes. Section III presents our optimal encoding/decoding 
framework for rate-diverse wireless network coding. Section 
IV describes the implementation of our encoding/decoding 
framework using LDLC. Section V proposes the BICM 
scheme for rate-diverse wireless network coding. Section VI 
provides simulation results. Section VII concludes this paper. 
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Fig. 1. The system model of the broadcast channel with side-information at 
the users. 

II. System Model 

We consider a network-coding assisted wireless broadcast 
problem. The system model is shown in Fig. [T] where we 
have a base station (BS) and two users (A and B). The BS 
wants to transmit different messages to users A and B. The 
message targeted for user A (B) is denoted by a vector of 
binary information bits S {0,1}^"^ (ms G {0,1}^^), 
where La and Lb the lengths of the vectors. User A (B) 
has side information ms (niyr), the message targeted for user 
B (A). This is a common scenario in the broadcast phase of 
the two-way relay channel based on physical-layer network 
coding Q; and this is also a special case of index coding 
systems Q. 

The BS can employ network coding to minimize the re¬ 
quired transmission time Besides network coding, the 

BS also needs to perform channel coding and modulation to 
generate the channel symbols that will be transmitted over the 
wireless channel. For simplicity, this paper considers a real¬ 
valued signal model. The extension to a complex-valued signal 
model is straightforward. We denote the vector of channel 
symbols transmitted by the BS by x^rc = / (rtiyt, ms), where 
xatc G consists of N channel symbols, and the function 
/ (•) incorporates the combined operation of channel coding, 
modulation and network coding. Thus, the data rate for user 
A (B) is Ra = La/N (Rb = Lb/N) bits per channel use. 
We impose an average power constraint Px on the channel 
symbols, i.e., E\\x.nc\\^ jN < Px- 

We model the wireless channel between the BS and user 
u G {A,B} as an additive white Gaussian noise (AWGN) 
channel with the path-loss effect: 

y„ =/SuXATc + n„ (1) 

where riu G is a vector of i.i.d. mean-zero, variance- 
Gaussian white noise components, and 0 < /3„ G K is 
the channel gain that models the path-loss effect between the 
BS and user u. The two channel gains 13a and Pb are likely 



Fig. 2. The capacity region of the broadcast channel with side-information 
at the two users. 

different due to the different distances of the users from the 
BS. 

After receiving ya (ys), user A (B) estimates its target 
message (ms) using ya (yb) and its side information 
(m^). Specihcally, we express the estimated target mes¬ 
sages as = gA{YA,mB), rhs = gBi.YB,mA), where 
gA{') and ps (•) denote the combined inverse operation of 
channel coding, modulation and network coding at users A and 
B, respectively. Henceforth, for brevity, we will simply call 
/ (•) the encoding scheme, and (•) and gB {■) the decoding 
schemes. Given the average power constraint over the channel 
symbols of BS, we ask the following two questions: 

i) What are the data-rate limits for the BS to reliably deliver 
messages to the two users? 

ii) What encoding/decoding schemes can be used to achieve 
these limits? 

We hrst identify the data-rate limits for the channel con¬ 
sidered here. If we just focus on the point-to-point single-user 
channel between the BS and one particular user u, the Shannon 
channel capacity, C„ = (1/2) log 2 (1 + SNRu), is the upper 
limit of the data rate for which reliable communication is 
possible as N ^ oo, where SNR^ = PxPujcr^ is the SNR at 
the receiver of user u. Considering the broadcast channel with 
side-information at the two users, references proved 

that as long as the date-rate pair {Ra,Rb) is within the 
capacity region given by {{Ra, Rb) ■ Ra < Ca, Rb < Cb}, 
the users can decode their target messages with arbitrarily 
small error probabilities. The capacity region is shown in 
Fig. 0 where we assume Pa > Pb- Obviously, the capacity 
pair {Ca, Cb) is the optimal data-rate pair that simultaneously 
maximizes the data rates for both users. 

This paper focuses on the second question: the encoding 
and decoding schemes to achieve the optimal data-rate pair 
(the capacity pair). Let us first consider the simple special case 
where the two channel gains are equal: Pa = Pb- Now, the 
point-to-point channel capacities of the two channels are equal: 
Ca = Cb = C, and the capacity pair becomes {C, C). For this 
rate-equal case, the encoding scheme for capacity achieving 
is rather straightforward. The first step is the linear network 
coding over the binary information: mjvc = niyi 0 ms, 
where 0 denotes the bit-wise XOR operation. Then, the 
network-coded message niArc is fed into a single-user channel 
encoder and modulator. At the receiver side of a user, a single- 
user decoding scheme can be used to obtain ihxc, the 
estimated target message is given by the bit-wise XOR of 
mxc and the side-information. As long as the data rate of 
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the used encoding scheme R can achieve the point-to-point 
single-user channel capacity C, the above simple encoding 
scheme can also achieve the capacity pair {C, C) for the rate- 
equal wireless network coding. 

Of interest to our paper here is the general rate-diverse 
wireless network coding where Pa ^ Pb and Ca ^ C's- 
How to achieve the capacity pair {Ca,Cb) for the rate- 
diverse case is not as obvious as the rate-equal case. In |[8), 
random coding is employed to derive the capacity region of 
the general probabilistic broadcast channel p{yA,yB Ixatc) 
with side information at the users. Then, considered the 
use of linear codes to achieve the optimal capacity pair for 
hnite-alphabet channels (channels with discrete outputs). By 
contrast, our paper here focuses on power-constraint AWGN 
channels. We investigate the encoding scheme and practical 
decoding scheme with affordable complexities to achieve (or 
to closely approach) the capacity pair. In particular, we put 
forth a design principle for rate-diverse wireless network 
coding. The design principle, referred to as the principle of 
virtual single-user channels, aims to transform the rate-diverse 
broadcast channel to two single-user channels states through 
encoding/decoding designs, thereby achieving the capacity 
pair which are basically the single-user capacities for user A 
and user B. 

III. Nested Lattice Codes Based Framework eor 
Achieving Capacity Pair 

This section describes an encoding/decoding framework for 
achieving the capacity pair of rate-diverse wireless network 
coding. The framework is based on nested lattice codes that 
have linear structures. We hrst give a preliminary on lattices 
and nested lattice codes in Section III.A. Then, in Section 
III.B, we show how nested lattice codes and its decoding can 
be used in our encoding/decoding framework to achieve the 
capacity pair. 

A. Preliminary on Lattices and Nested Lattice Codes 

A real lattice A of dimension iT is a discrete subgroup 
of (K < N) closed under addition and reflection: if 
Ai,A2 S a, then Ai -f A 2 G A; if A G A then —A G A. 
The lattice points (vectors) of A are generated by taking all 
integer linear combinations of K independent basis vectors. 
The K basis vectors can be written into an N x K generator 
matrix G G Therefore, a lattice A is specified 

by its generator matrix G and can be always written as 
A(G) = {A = Gb:bGZ^}. 

Given a lattice A and a lattice point A G A, the Voronoi 
region of A is defined to be the set of all vectors in 
that are closest to the lattice point A: V (A, A) = 
{x G : ||A — x|| < ||A' — x|| , VA' G A, A' ^ A}. The 
Voronoi region of A = 0 is called the fundamental 
Voronoi region of the lattice and it is denoted by V. 
Since every vector in R^ can be uniquely written as 
X = A -I- r, where A G A and r G V, a lattice quantization 
operation is a function that maps a vector x G R^ to a 
lattice point of A according to the minimum Euclidean 
distance rule: Qa (x) = argmin^'eA ||x — A'|| = A; a 


lattice modulo operation is to get the quantization error: 
X modA = X — Qa (x) = r. 

The volume of V is denoted by Vol (V) and it can be shown 
that Vol (V) = fydx= i/det (G^G). The second-order mo¬ 
ment of V is defined as (V) = ;^£’||U||^ — W fv wfv) 
where U is a random vector uniformly distributed over V. 
The normalized second-order moment of V is defined as 
G (V) = ^ which is invariant to the scale of the lattice 

(i.e., if the generator matrix G were to be scaled by a non-zero 
scalar, G (A) would remain the same). The above definitions 
of volume and (normalized) second-order moment are not only 
restricted to Voronoi regions of lattices; they are applicable to 
any regions over the space of R^. In other words, in place of 
V, for any subset of S' C R^, we will also make use of the 
generalized definitions for the second moment cr^(S) and the 
normalized second moment G(S) later in this paper. 

If the lattice A' is a subset of another lattice A, A' C A, 
we say A' is nested in A. A pair of lattices (A', A) is called 
a nested pair if A' C A, where A' is called the coarse lattice 
and A is called the fine lattice. For example, (gZ^,Z^) is a 
nested pair, where g is a non-zero integer. 

We now apply lattices to the coding problem. In this 
subsection, we consider a point-to-point single-user channel: 
y = /3x -f n, where x is the vector of the channel symbols, 
y is the vector of the received signals, /3 is the channel gain, 
and n is the vector of i.i.d. real AWGN components with 
mean-zero and variance-cr^. All vectors here are length-A^ 
vectors. The vector of the channel symbols x is the codeword 
for conveying information message m and it is subject to 
an average power constraint E\\x\\^ jN < Px- The aim is 
to achieve the channel capacity G = (l/2)loen (1 -I- SNR), 
where SNR ^ Px0^ jal- 

For lattice coding, the codeword x is a point chosen from 
a specific lattice A that contains an infinite number of lattice 
points. We call A the coding lattice. The power constraint 
on X means that only a finite number of the lattice points 
of A can be chosen as the codewords. To limit the power of 
codewords, we can choose a subset of lattice points within 
a certain region to be the codewords. This region, typically 
centered on the zero lattice point, is referred to as the shaping 
region. Overall, we need to take into account the following 
two aspects for designing lattice codes Q. 

• The granularity of the coding lattice A is represented by 
its fundamental Voronoi region V. The volume Vol (V) 
determines the inter-codeword Euclidean distance, hence, 
it determines the decoding error probability. 

• The structure of the shaping region determines the power- 
volume tradeoff, hence, the gap from the channel capac¬ 
ity. 

Two key questions are what is a good lattice for coding 
and what is a proper shaping region that satisfies the power 
constraint. 

To see what is a good lattice for coding, let us first remove 
the power constraint on the codeword x. In this case, since 
the transmission power as well as the data rate is infinite, 
any point of a lattice can be chosen as the codeword. At 
the receiver side, the maximum likelihood (ML) decoding is 





OPTIMAL RATE-DIVERSE WIRELESS NETWORK CODING 


5 


employed to search for a lattice point nearest to the received 
vector. Obviously, the decision regions of the ML decoding 
are the Voronoi regions and this ML decoder essentially is 
the lattice quantizer Qa ( )■ The performance of the code is 
expressed by the decoding error probability Pr (rfi 7^ m). 
Since decoding errors occur when the noise vector goes 
beyond the Voronoi region of the transmitted lattice point, 
the decoding error probability is Pr (m 7^ m) = Pr (n ^ V). 
The decoding error probability is determined by the “shape” 
of V, which depends on the lattice being used, and the 
size of V relative to the noise power, which can be 
expressed as the volume-to-noise ratio (VNR) given by 
7 (A, cr^) = Vol(V)^^'^y/tr^. A good lattice code should 

have low Pr (n ^ V) for a fixed 7 (A, cr^). According to |j^, 
we have the following result for lattices good for coding. 

Goodness of Lattices for Coding: A sequence of lattices 
indexed by their dimension is said to be good for coding if 

i) for a target decoding error probability Pr (m f m), 
where 0 < Pr (in f m) < 1, VNR 7 (A, cr^) required 
to achieve the target Pr (in f m) approaches 27re as N 
goes to infinity (i. e., limAr_>.oo 7 (A, cr^) = 27re); and 

ii) for a fixed 7 (A, cr^), where 7 > 27re, the decoding error 
probability Pr (in f m) vanishes exponentially in N 
that is greater than 27re, Pr (in f m) vanishes exponen¬ 
tially in N (i.e., limAr->.oo Pr (in m) = 0). ) 

Reference | |22l showed such lattices exist. 

We next examine the power-constraint case to see how 
to choose a good shaping region. The (normalized) second- 
order moment is a metric for the average power of a random 
vector uniformly distributed over a given shaping region, thus 
it measures how good the shaping region is. Among all N- 
dimensional bodies of a fixed volume, the body with the mini¬ 
mum normalized second-moment is the N-dimensional sphere, 
denoted by The normalized second-order moment of 

the W-dimensional sphere decreases monotonically 

with N and approaches the limit l/27re as N ^ 00 0 , 0 . 
Thus, the optimal choice for the shapin g re gion is the N- 
dimensional sphere denoted by 1 ^, 

To introduce the notion of shaping loss, we consider a 
simple choice for the shaping region: the A^-dimensional hy¬ 
percube. It is well-known that there is a shaping loss between 
the average powers of hypercube shaping and sphere shaping 
10. The hypercube is the fundamental Voronoi region of 
the lattice Z^, and its normalized second-order moment is 
G (^) = 1/12, \/N. Therefore, compared with the optimal 
sphere shaping, the shaping loss for the hypercube shaping 
is 7 s (Z-^) = (l/27re)/G (Z^) = 7re/6 as N ^ 00 . This 
7 re/6 (1.53 dB) shaping loss for the hypercube shaping dose 
not vanish as the SNR increases, and thus it is an undesired 
feature, especially in the high SNR regime |j0. 

To maintain the optimalities for both coding and shaping, 
| [24] showed that a spherical lattice code (the intersection of a 
lattice good for coding with an A^-dimensional sphere of radius 
V NP) can arbitrarily approach the channel capacity. Although 
the sphere region can maintain the shaping optimality, such 



Fig. 3. A nested lattice code with N = 2 and Vol (Vs) = I6V0I (V) (the 
shaded region is the used codebook). 

spherical lattice code destroys the linear structure of the orig¬ 
inal coding lattice A. Moreover, the optimal ML decoding for 
such spherical lattice code is not the lattice quantization Qa (•) 
anymore, because the ML decoding regions for codewords 
are not identical and some are not bounded. By contrast, in 
the unconstraint case, the lattice quantization Qa (•) used for 
ML decoding ignores the boundary of the code and preserves 
the symmetry of the lattice structure in the decoding process: 
it is much less complex as far as decoding is concerned. 
For practical implementation, it is desirable to preserve the 
linear structure of lattices both in the encoding and decoding 
processes. 

Ref. developed a lattice framework that can reliably 
transmit at rates up to the channel capacity. This framework 
is called nested lattice codes and its general idea is to 
make use of a nested pair of lattices (As,A), where the 
coarse lattice A^ is used for shaping and the fine lattice 
A is used for coding. We denote the fundamental Voronoi 
region of the coarse lattice Ag by "Fg, and the volume of Vg 
by Vol/Fg). The nested lattice code is generated by taking 
the intersection of the fine lattice used for coding with the 
fundamental Voronoi region of the coarse lattice used for 
shaping: C = {A n )2g}. The coding rate of the nested lattice 
code is i? = jTiog 2 \C\ = ^loga Fig. illustrates an 

example for the codebook of a nested lattice code, where the 
dimension is = 2 and Vol ("Fg) = I 6 V 0 I /F) (thus, R — 2). 
We state the following result for the goodness of lattices for 
shaping. 

Goodness of Lattices for 5/iapin^ A sequence of lattice 
aI'^^ is good for shaping if lim G {a^s^^\ = l/27re. Such 
lattices exist as shown in | |25| . 

It is known that the normalized second-order moment 
of a lattice is always larger than l/27re, the normalized 
second-order moment of a sphere with infinite dimensions. 
The goodness of lattices for shaping indicates that as the 
dimensions become sufficiently large, there are lattices whose 

' In lattice literature, this feature is also termed as the goodness of lattices 
for MMSE quantization. 
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fundamental Voronoi region approach a sphere as their 
normalized second-order moments go to l/27re. Therefore, 
such lattice is asymptotically optimal for shaping in terms 
of that its normalized second-order moment approaches 
l/27re. The authors of pT) showed that nested pair of lattices 
(As, A), where the coarse lattice As is good for shaping 
and the fine lattice A is good for coding, exist for any 
required rate. Therefore, based on such nested pair of lattices 
(As, A), the channel capacity can be potentially be achieved 
in all SNR regimes. However, before that, there are still two 
important ingredients of nested lattice codes; the dithering 
operation and minimum mean square error (MMSE) scaling 
at the encoding and decoding processes. We now give the 
complete description for the encoding and decoding processes 
of nested lattice codes ID- 

Encoding and Decoding Processes of Nested Lattice Codes: 

• Encoding: First, the message m is mapped to a codeword 
c = (p (m), where f (•) is the message-to-codeword map¬ 
ping function, the codeword x belongs to the codebook 
of the used nested lattice code C = {A n Vs} . Then, the 
transmitted vector is generated according to 

X = [c — d] mod As (2) 

where d G Vs is the dithering vector that is uniformly 
distributed over the shaping region Vs- The dithering 
vector d is known at both of the encoding and decoding 
processes. 

• Decoding: The estimate for the transmitted codeword is 
computed according to 

c = QA([ay' + d])mod As (3) 

where y' = /3“^y is the channel-gain-normalized 

received vector, a = SNR/{1 + SNR) is the MMSE 
coefficient used to scale the channel-gain-normalized 
received vector y' before sending it to the lattice 
quantizer for decoding. 

The above encoding and decoding processes of nested lattice 
codes are illustrated in Fig. The dithering operation in 
© can ensure that the distribution of x is the same as that 
of d (c.f. Lemma 1 of Q). Therefore, as long as the used 
shaping lattice is scaled to have the second-order moment 
Px, the power of the transmitted vector of symbols x is Px- 
Furthermore, if the shaping region Vs approaches a sphere 
as iV —> oo, X will have white Gaussian distributions, as 
desired by Shannon theory. The use of MMSE scaling in 
Q plays an important role for the purpose of achieving the 


channel capacity, especially in the low SNR regime. Please 
see Q for more details about the MMSE scaling. In Q, it 
was proved that the lattice quantizer decoding in Q suffices 
to be optimal. We end this preliminary on nested lattice codes 
here. The reader is referred to ||9) for further details. In 
conclusion, with nested lattices codes, the channel capacity 
C = (1/2) log2 (1 + S'A^i?) can be achieved in all SNR 
regimes. In the next section, employing nested lattice codes, 
we will present a framework for achieving the capacity pair 
of the rate-diverse wireless network coding. 

B. Nested Lattice Framework for Achieving Capacity Pair 

Consider the rate-diverse network coding problem. The two 
users have different channel qualities, thus different channel 
capacities. As a consequence, the codes operating at different 
channels have different rates. To develop two nested lattice 
codes with different rates, we employ two nested pairs of 
lattices {As,Aa), {As,Ab), where the two different fine 
lattices A^i and Ab are used for the coding of user A 
and B respectively, the same coarse lattice A^ is used for 
the shaping. The corresponding two nested lattice codes are 
Ca = {Aa n Vs}, Cb = {A_b n Vs}. We employ Ca, Cb in 
our nested lattice framework to achieve the capacity pair of 
the rate-diverse network coding problem. 

The proposed encoding scheme /(•) at the transmitter of 
BS first maps the messages m^, ms into the codewords: 
CA = (pA (m^i), Cb = 4>b (ms), where c^i G Ca, cb G 
Cb, and (•)’ 4>b {■) are the message-to-codeword mapping 
functions for codes Ca, Cb- Then, we perform the network 
coding operation over the codewords to form the network- 
coded codeword 

cnc = [ca + Cb] mod A^ (4) 

Finally, like nested lattice codes for point-to-point channels, 
we perform dithering operation to generate the vector of 
channel symbols 

xatc = [cnc - d] mod A^ (5) 

Since cxc G Vs, the distribution of xjvc is still a uniform 
distribution over the shaping region As (the same as the 
dithering vector d). Thus, if As approaches a sphere as the 
N grows, the transmitted vector xxc will look like a white 
Gaussian noise. This satisfies the requirement on the channel 
inputs by Shannon theory. Using the distributive law of the 
modulo arithmetic, the operations © and © can be combined 
into 

xatc = [[ca + Cb] mod As - d] mod As 
= [ca -I- Cb - d] mod As 




























OPTIMAL RATE-DIVERSE WIRELESS NETWORK CODING 


7 



Based on the expression in we illustrate the proposed 
encoding scheme for rate-diverse wireless network coding in 
Fig. ID 

The decoding process at the receiver of user A, qa {•), 
performs the following steps in sequence; the network decod¬ 
ing (the substruction of its side-information), the scaling and 
de-dithering, and the lattice quantization with respect to the 
coding lattice A^. Thus, the estimate for the target message 
of user A is given by 


/ \ 



/ 

p p 


Ca = Qaa 


<aa {Pa^Ya - cb) +d 




Step I iNetwork Decoding 



Step II rScaling and Dedithering 


Step III :Lattice Quantization 

(7) 

where 01,4 = SNRa/{^ + SNRa) is the MMSE coefficient 
for user A. User B performs a similar decoding process. The 
decoding processes for rate-diverse wireless network coding 
are also illustrated in Fig. 

Substituting the received vector at the user A, = 
Pa^nc + Ha, into 0 and making some manipulations give 


ca = Qka ([“aCa - o-aA + aAl^A^^A + d]) modAs 
= Qa^ {[a AC a - a^d -f uaPa^taa + d] modA^ 
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( 

\ 


\ 

Qaa 


OLA 

[cA 

- d] mod As +P2^nA 

+ d 

mod As 





A 







-V- ^ 




\ 


1 

^ J 


/ 


= Qaa {[o^aY'a + d ] modAs) 

= Qaa {{<AAy'A + d]) modAs 

( 8 ) 

where xa = [ca — d] mod Ag that is the vector of virtual 
point-to-point channel symbols obtained by performing single- 
user encoding on message niA using nested lattice code Ca 
designed for user A, and y'^ = xa + /3)^^nA is the equiv¬ 
alent channel model for user A after the network decoding 
operation. 


In (j^, we can see that the equivalent channel left for user 

A, y'A — XA + /S^^riA, is the same as a point-to-point 
single-user channel where message niA is conveyed by vector 
Xa. We also note that the decoding operation for rate-diverse 
wireless network coding expressed by (j^ has the same form as 
the decoding operation in a point-to-point single-user channel 
expressed by 0 - Compared to the point-to-point single-user 
channel from BS to user A, the SNR of the equivalent channel 
is still SNRa = Px^a/’^u (Fe. the SNR is not reduced in 
the rate-diverse wireless network coding case). In this sense, 
user B is totally transparent to user A. Therefore, if the used 
nested pair (As,Aa) is good for shaping and coding for the 
point-to-point single-user channel of user A, we can achieve 
the channel capacity Ca using nested lattice code Ca- With 
the same decoding scheme, the same result holds for user 

B. Therefore, as long as nested pairs (As,Aa), {As,Ab) 
are both good for shaping and coding, the above proposed 
encoding/decoding framework can achieve the capacity pair 
(CajCb). Such nested pairs do exist and can be obtained 
using the construction method proposed in | [26| . 

With the linear structured nested lattice codes, the above 
encoding/decoding framework is optimal for achieving the ca¬ 
pacity pair of the rate-diverse network coding. The framework 
brings out a general design principle for the encoding/decoding 
processes. The principle, referred to as the principle of virtual 
single-user channels, can be stated as follows: 

Principle of Virtual Single-User Channels ; 

i) At the transmitter side, the encoding scheme performs 
coding for the two users separately, using the codes 
with different rates that matches with the capacities of 
their respective channels. The two codewords are then 
combined by a network coding operation. 

ii) At the receiver side, each user applies network decoding 
on the received signal so as to extract the signal contain¬ 
ing only the desired codeword. The signal of the other 
user should be totally transparent to the user after the 
network decoding process. In particular, the remaining 
signal is exactly the same as that of a single-user point- 
to-point communication, including the effective signal 
strength and the effective noise thereof. Applying single- 
user decoding to extract the desired message has perfor¬ 
mance exactly the same as that in conventional point-to- 
point communication. 
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The above nested-lattice-code framework adheres to this 
design principle. As will be seen later, this principle can 
also be realized using other codes with simpler decoding 
algorithms. 

Although the nested-lattice-code framework is optimal the¬ 
oretically, its exact implementation faces many difficulties. In 
particular, the lattice quantization decoding, which searches 
the transmitted lattice point over the lattice space, has unaf¬ 
fordable complexity as the codeword length N grows. For the 
existing lattice quantization decoding methods, the complexity 
increases exponentially with N. There have been efforts to 
design practical lattice codes with implementable decoding 
algorithms p7|. I n the next section, we employ one of 
such codes, LDLC pO), to realize the design principle. In the 
section after that, we show that the same design principle can 
also be realized by BICM. 

IV. Implementation oe the Proposed Framework 
Using LDLC 

We provide the background on LDLC in Section IV.A. 
We apply LDLC in our rate-diverse wireless network coding 
framework in Section IV.B. 


A. Low Density Lattice Codes (LDLC) 

As discussed in Section III, lattice codes use the points of a 
specific lattice as codewords (messages). The codewords of a 
lattice code are generated by x = Gb, where x G is the 
codeword, G S is the generating matri)(|^ and b S 

is the integer message vector. LDLC is a type of lattice codes 
where the parity check matrix, defined as the inverse of the 
generating matrix H = G“^, is sparse | [l^ . 

The number of non-zero elements in the i-th row (column) 
of H is called the degree of the i-th row (column). The sparse 
structure of H implies that the degrees of its all rows and 
columns are far less than N. For LDLC, the non-zero elements 
of H can adopt values in K, in contrast to ordinary LDPC 
where the elements of H are typically elements of a finite 
field. 

In practice, the elements of the integer message vector b 
only can be selected from some finite constellations, rather 
than the infinite Z^. We denote the i-th element of the integer 
message vector h hy bi G Ci, where Ci is a subset of integers. 
The coding rate of such LDLC is R = \^i\^ 

where \Ci\ is the cardinality of Ci. 

In pO) , Latin square LDLC were considered. Latin square 
LDLC are defined by imposing some restrictions to their 
parity check matrixes. For a Latin square LDLC, all rows and 
columns of H have the same degree d; a sorted sequence 
of d non-zero values, {hd > hd-i > ■ ■ ■ > hi > 0}, is the 
generating sequence of H; the d non-zero elements in one 
row or column of H are {hd, hd-i, ■ ■ ■, hi] with a change 
of positions and signs for different rows and columns. Given 
a generating sequence, the parity check matrix of a LDLC 

^We assume square LDLC generating matrices in this section. The results 
can be easily extended to non-square generating matrices. 


can be constructed using the proposed method in pO) . For 
example, {1,0.8, 0.5} is a generating sequence and 
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is the corresponding parity matrix of a LDLC with degree 
d = 3 and dimension N = 6. 

For communications over power-constrained AWGN chan¬ 
nels, only codewords within some shaping region can be cho¬ 
sen as the valid codewords. Directly transmitting codewords of 
LDLC without shaping will increase the average transmission 
power. Fig. (a) illustrates this situation using an example. 
In Fig (a), we generate a two-dimensional lattice using the 

generating matrix G = ^ ^ v^/2 ) ' black dots are the 

lattice points. The squares are the lattice codewords selected 
by X = Gb, where Ci = {—2, —1,0,1, 2}for all i. We have 
25 lattice points are selected as the codewords, and the coding 
rate thus is i? = log225/2 = 2.3. Fig. (a) also shows 
the nested lattice shaping region. From Fig. (a), we can 
intuitively see that the codewords outside the shaping region 
will greatly increase the average transmission power. 

To reduce the average transmission power, we must map all 
the codewords x = Gb into the codewords x' that are within 
the shaping region. This shaping operation can be expressed 
by mapping integer message vectors b to some other integer 
vectors b', such that the obtained codewords x' = Gb' lie 
within the shaping region. Note that b/ may not be an element 
of Ci anymore. We will describe how to perform shaping for 
LDLC in next subsection. 

Due to the sparse structure of H, a belief propagation 
(BP) algorithm can be used to decode the LDLC codeword 
at the receiver side. Given the parity check matrix H, we can 
construct the factor graph that models the constraint imposed 
by LDLC encoding. Then, we derive the BP decoding as a 
message passing algorithm on the factor graph using the sum- 
product rule p^ . Since the factor graph for LDLC contains 
loops, the BP algorithm will need to go through several 
iterations. Although the BP decoding algorithm is in general 
not optimal for lattice decoding, its linear complexity is its 
advantage over more complicated decoding in practice. Com¬ 
pared with BP decoding for LDPC, BP decoding for LDLC 
computes continuous-valued messages rather than discrete¬ 
valued message; the check nodes of the factor graph force the 
summation of their incoming messages to integers with certain 
pre-defined range. For details on BP decoding for LDLC, we 
refer the readers to ID- 


B. Wireless Network Coding Using LDLC 

We now apply LDLC to the rate-diverse wireless network 
coding problem. The codewords ca, cb transmitted to user 
A and B are generated by applying LDLC encoding: ca = 
GyibA, Cb = Gsbs, where b^ is the integer message vector 
of user u, G„ is the LDLC generating matrix of user u. After 
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(a) (b) (c) 


Fig. 6. An example of the coding and shaping aspects for single-user LDLC and wireless network coding LDLC: (a) Single-user LDLC with coding rate 
2.32; (b) Single-user LDLC with coding rate 1.59; (c) Wireless network coding LDLC obtained by directly adding together the lattice codewords of user A 
in (a) and that of user B in (b). 


network coding and lattice shaping, the transmitted codeword 
is Ctj. Note that the dithering vector d in (6) usually is not 
used in real systems, and thus we do not consider it hereafter. 

The lattices A a and Ab generated by and Gb are 
the coding lattice for user A and B, respectively. For shaping 
purposes, our optimal framework requires a coarse lattice 
that is simultaneously nested in the two coding lattices, i.e.. 
As C A^, As C Ab. The code rate for user A (B) is 
determined by the intersection between A a (Ab) and Vs- The 
inverse matrices of G^ and Gb must be sparse to allow 
practical BP decoding. The generating matrices G^, Gb, Gs 
of such A^, Ab, As can be found as follows. 

We hrst hnd an LDLC generating matrix G using the 
construction method proposed in pO) . The inversion of G 
is sparse. We call G the basic generating matrix. Different 
coding rates are assigned to the two users by applying different 
constellations to the integer message vectors of the two users, 
i.e., bu,i G Cu,i = {-Lu,i,-Lu,^ + l,■■■ where 

integer ^ is the i-th element of b„, ^ is the constellation 

for bu,i, and ^ is an integer parameter that determines the 
corresponding constellation size via \Cu,i \ = 2Lu,i. Note that, 
for the same user u, we may assign different constellations 
to different integer message element ^ To limit the av¬ 
erage transmission power, we choose the shaping lattice as 
Gs = GM, where M = diag {Mi, M 2 , ■'' ^ is an 
N X N diagonal matrix. The z-th main diagonal element 
of M, Mi, is set to be the least common multiple (LCM) 
of 2LA,^ and 2LB,i: M, = LCM {2LA,^,2LB,i}■ Finally, 
we choose the coding lattice for user u as G^ = GM„, 
where M,i = diag {Mua, AIu, 2 , ■ ■ ■, AIu,n} is an N x N 
diagonal integer matrix whose i-th main diagonal element 
is set as i With the above setup, we can 

easily check that As C Aa, As C Ab are satished — the 
shaping optimality of the nested lattice framework for rate- 
diverse wireless network coding is guaranteed. 

Fig. 0 (a) and Fig. (b) show examples for LDLC lattices 
with two dimensions (N — 2) for users A and B, respectively. 
The black dots are the lattice points generated by the basic 
generating matrix G. The squares in Fig. (a) are the 
lattice codewords for user A selected by applying constellation 
^A,i = {—2, —1, 0,1, 2} to bA,i for all i. The squares in Fig.|^ 


(b) are the lattice codewords for user B selected by applying 
constellation Cb,* = {~1,0; 1} to 6 b, i for all i. The code rates 
for users A and B are 2.32 and 1.59 respectively. In Fig. 

(c) , the red squares indicate the lattice codewords obtained by 
directly adding together the lattice codewords of user A in 
Fig. 1^ (a) and that of user B in Fig. (b) without shaping. 
From Fig. (c), we can see that, without shaping, the average 
transmission power increases quite heavily compared with the 
single-user case, due to the addition. Therefore, shaping is 
more crucial for the network coding case. The next question 
is how to perform the nested lattice shaping for rate-diverse 
wireless network coding using the coarse lattice A^. 

As long as the coarse lattice Ag (equivalently, Gg) used 
for shaping is specified, the lattice shaping operation, = 
[c^ + cb] modAg, can be expressed as 


^NC = ca + Cb — Qa, (ca + Cb) 
= ca -f Cb - Ggk 


(9) 


where Ggk = Qa^ (ca -t- cb), k G , is the coarse lattice 
point that is closest to ca -f cb. Substituting ca = GAbA = 
GMAbA, Ca = GBbB = GMBbB and Gg = GM into 
(0, we have 


xtvc = GMAbA -f GMBbB — GMk 

/ \ 


= G 


M^AbA -f M^BbB —Mlk 


A. 

= OnC 


= b'jvc 


NC 


( 10 ) 


where h 'nc = i^NC — Mk with hjAc = MAbA -t- MBbB. 
The shaping operation in ( [T0) i is to map the network-coded 
integer message vector b^vc to another network-coded integer 
message vector h' nc = b^vc ~ Mk, such that the resulting 
coding lattice point (the codeword xjvc = GbVc) lies within 
the fundamental Voronoi region of the shaping lattice. 

The target of the shaping operation is to solve the quantiza¬ 
tion problem Ggk = Qa^ (ca + cb), which is equivalent to 
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finding the integer vector 

k = argmin ||cA + c_B - Gsk'll^ ( 11 ) 

k' 

Since the inverse of = GM is still sparse, at first sight, 
we can also employ the BP algorithm developed for decoding 
LDLC to solve the shaping problem 0 approximately. How¬ 
ever, |T§ showed that the BP algorithm cannot give satisfac¬ 
tory shaping performance. The reason is that the distribution of 
the effective noise (i.e., the deviation from G^k) in the shaping 
problem (i.e., the codeword c .4 -f cs) is a uniform over the 
space of codewords, not a Gaussian distribution such as that of 
channel noise. Due to the uniform distribution, the probability 
of the codewords ca + cb appearing at the boundaries of the 
Voronoi regions of the shaping lattice will not decrease with 
the distance between the codewords and the centers of the 
Voronoi regions. Therefore, the variance of the effective noise 
in the shaping problem is much larger than that in the decoding 
problem. The BP algorithm that iteratively and approximately 
solves the decoding problem fails to reduce the power of xjvc' 
as a result. 

Alternatively, | |28t employed the M-algorithm to solve the 
shaping problem for LDLC in point-to-point channels. For the 
shaping problem of LDLC in rate-diverse wireless network 
coding, we also use the M-algorithm. 

The M-algorithm is a sequential tree search algorithm 1^ . 
To apply it to the shaping problem ( [TT] i, we first perform QR 
decomposition on the LDLC parity check matrix = Q^R, 
where is an orthogonal matrix and R is a upper triangle 
matrix. Then, H can be written as H = TQ, where T = R^ 
is a lower triangle matrix. Multiplying H = TQ on both 
sides of Xatc = Gh'^vc^ we have the following new check 
relationship; 

TxAfc = bVc = bwc — Mk (12) 

where Scjsic = Q^nc is the codeword after the transformation 
by Q. Since T is a lower triangle matrix, the check relation¬ 
ship ( [T 2 | ) can be decomposed into N check equations 

i— 1 

TijXNCJ + Ti^iXNC,i = bNC,i — M^ki, 

j=i ttj; 


where Tij is the {i,j)-th element of T. The i-th check 
equation in ([T^ corresponds to the i-th row of ( [T^ . Originally, 
the shaping problem in (11 1 aims to minimize ||xArc||^- Since 
xatc = Qxatc and Q is an orthogonal matrix that will not 
change the power of knc, we now change our objective 
to minimizing |jx 7 vc||^ to solve the shaping problem. The 
triangular structure of T in (12i suggests that we can make 
use of a tree search method, such as the M-algorithm, over all 
possible vector k to find the one that minimize the power of 
XatC- 


We consider a tree rooted at a dummy node. The nodes at 
the i-th level of the tree are labeled by different sequences 
of ki,k 2 , - ■ ■ ,ki- We note that given a particular sequence 
ki,k 2 , - ■ ■ ,ki, the sequence xnc,i,xnc, 2 , ■ ‘ ‘ ,xnc,i is also 
fixed due to ([T3l). Each node ki,k 2 , - ■ ■ is associated with 
a metric W^ncj li^- The optimal ML tree search requires 


tracing all the paths on the tree. The number of possible paths 
on the tree increases exponentially with the depth of the level. 
Since the depth of the tree, the length of k, is large for 
non-trivial lattice codes, this may incur huge computational 
complexities. Therefore, sub-optimal algorithms with lower 
complexity, but with good performance, are needed to perform 
the tree search task for LDLC shaping. The sub-optimal M- 
algorithm starts the tree search at the root node, and visits 
the levels from the one to the next one. At every level of the 
tree, the M-algorithm just retains M paths rather than all the 
possible paths, where M is a design parameter. The node at 
the i-th level ki,k 2 , - ■ ■, ki-i, ki is a child node of its parent 
node fci, fc 2 , • • •, ki-i at the {i — l)-th level, and the metric of 
the node ki,k 2 , - ■ ■, fci-i, ki is recursively computed as 


\\xNC,j\f 

i=i 


i-1 


\\xNC,j\f + \\xNC,i\? 


where 


XNC,i 



^ ^iki I 

i=l / 


(14) 


(15) 


is the i-th element of xjvc computed according to We 
assume h S {- \Mi/2\ , - \Mi/2\ -f 1, • • •, [Mi/2J - 1}. 
Among all the paths reached at the z-th level, we just retain 
the M paths that have the least M values of the metric 
l|2?AfC,jiy Then, the search continues to the (i-f 1)- 
th level of the tree. After the search finishes visiting the leaf 
nodes of the tree, we obtain M paths labeled by the sequence 
ki,k 2 , -' ‘ i^N, and finally we select the one with the least 
metric value ll^tvGjll^ as the final estimate for k. 

At the receiver of user A, a network decoding operation 
first substracts the side information from the received signal: 
y'A = Ya — I^aG'M.b^b- Then, BP channel decoding is used 
to estimate — Mk. After that, the estimate for 

the target message is immediately obtained by the modulo op¬ 
eration b^ = (b'^mod M), which can be performed in 
an element-by-element manner: bA,i = (6'yi_imod Mi) 
for all i. We remark here that the use of the modulo operation 
for recovering hA (b^i) avoids the need to derive k directly 
at the receivers. 

We remark here that the original constellation of ^ is 
^u,i = {-Lu,i, -Lu,i + 1, • • •, Lu,i - 1} ; after shaping, the 
constellation of b'^i is changed to Cu,i + Miki. For the 
decoding of the unshaped LDLC, the receiver just consid¬ 
ers the possible values of bu^ over Cu,i- For the decod¬ 
ing of the shaped LDLC, the receiver needs to know the 
constellation of b'^ j. Here, we can empirically choose a 
constellation for b'^ j. Typically, for the shaping operation, 
the possible values of ki is bounded by \ki\ < £ for all 
i, where e is a small positive integer. For example, in our 
simulations of Section VI, we observed that e = 2. There¬ 
fore, we can expand the constellation of 6^ j to ^ = 
{-Lua - sMi, - eMi -f 1, • • •, -f eMi - 1} and 

use it in the decoding]^ 


^ In practice, the transmitter of BS can embed e in the packet being sent 
as metadata to facilitate the decoding by receiver. 
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Since the principle of virtual single-user channels is sat¬ 
isfied, we predict that the performance of the above LDLC 
based rate-diverse wireless network coding scheme is the same 
as that of a LDLC based single-user point-to-point scheme. 
Moreover, as will be elaborated in the next section, there is 
an additional shaping gain for the user with the lower data 
rate. 

C. Additional Shaping Gain for Rate-Diverse Wireless Net¬ 
work Coding Using LDLC 

In the optimal lattice coding framework, the shaping opera¬ 
tion employs a dithering vector d to disseminate the transmis¬ 
sion power over the whole shaping region, as indicated by Q 
for the point-to-point case and by (j^ for the network coding 
case. The distribution of the dithering vector is a continuous 
uniform distribution. The distribution of lattice points is a 
discrete uniform distribution. The average power of the lattice 
points within a shaping region can be changed by varying 
the density of these lattice points. By increasing the density 
of the lattice points within the shaping region (piling more 
lattice points into the shaping region), the average power can 
approach its lower bound — the average power of a continuous 
and uniform random variable over the shaping region. After the 
shaping operation as in Q, the distribution of the codeword x 
is the same as the dithering vector d (a continuous and uniform 
random variable). Therefore, from a theoretical perspective, 
this dithering operation can reduce the transmission power of 
a nested lattice code. However, in practice, this dithering vector 
is not used in the coding process of real systems, since the 
generation and storage of the A^-dimensional random vector 
d G will impose a large complexity on the system as N 
grows. 

Without dithering, we perform the shaping operation for 
network coding as x^rc = [cA + csJmod Ag. Due to the 
network coding operation, the lattice points of user A (B) 
behave somewhat like the dithering vector for user B (A). 
The density of the codewords, Xjvc> will be the higher one 
between the densities of ca and c b ■ Therefore, compared with 
the point-to-point case, the user with a lower data rate (lower 
density of its lattice points within the shaping region) achieves 
an additional shaping gain. We can demonstrate this additional 
shaping gain mathematically using a simple one-dimensional 
lattice example. 

Consider two one-dimensional lattices; = Z for 

user A and Ab = 2Z for user B. The shaping re¬ 
gion is taken as [—L, L], where we assume L is a po¬ 
sition even number for simplicity. As a consequence, the 
codewords of the nested lattice codes for users A and 
B are ca G {—L, —L-f 1, • • •, L — 1, L} and cb G 
{—L, —L -\- 2, ■ ■ ■, L — 2, L} respectively. The average pow¬ 
ers of the codes for user A and B are Pa — E (c^) = 
L(L-fl)/3 and Pb = -E (c|) = L(L + 2)/3. After 
the operations of network coding and shaping; xnc = 
[cA + csJmod 2L, the transmitted codewords are xmc G 
{—L, —L -f 1, • • •, L — 1, L} and the average power of X]sic 
is Pnc = E {x^Aifj) = Pa = L{L -\- l)/3. This means that, 
compared with the point-to-point case, there is an additional 


shaping gain, lOlog^^PB/ Pnc = 101ogio(i + 2 )/(L-f 1), 
for user B, which has a lower data rate in the rate-diverse 
wireless network coding case. We can observe this additional 
shaping gain from the numerical results in Section VI. 

Recall that our design principle of virtual single-user chan¬ 
nels states that we should strive for a design which achieves 
the single-user performance. The implication is that single- 
user performance is some sort of an upper bound to be aimed 
for, and that the two-user rate-diverse system cannot have 
performance better than the best of single-user systems. The 
above result, at first sight, may appear to have contradicted the 
above statement. A more careful examination can resolve this 
“paradox” easily. Note that the additional shaping gain in the 
two-user rate-diverse system is caused by the denser effective 
codebook, which induces some sort of a dithering effect on an 
un-dithered system. If dithering had been used, then the denser 
codebook would have no effect. In other words, it would not 
possible for the two-user rate-diverse system to beat the single- 
user system if dithering had been applied in the latter in the 
first place (i.e., it is not possible to beat the best single-user 
system). 

V. Implementation oe the Proposed Framework 
Using BICM 

BICM is a framework for combining binary channel codes 
with high-order modulations (e. g., QAM) in- For the single- 
user point-to-point problem, BICM first encodes the binary 
information bits using a binary code, then interleaves the 
coded bits, and then maps the interleaved bits to modulated 
signals.. By varying the code rate of the binary channel code 
and the size of the modulation signal set, BICM can support 
various target data rates. Moreover, it has been shown that 
the achievable rate of BICM is very close to the single-user 
channel capacity HI, m- 

We next show how to apply BICM to the rate-diverse 
wireless network coding problem. We design the BICM system 
to adhere to the principle of virtual single-user channels 
expounded in Section III.B. Fig. shows the block diagram 
of our proposed BICM scheme. The BS first encodes the 
binary information bits of users A and B (m^ G {0,1}^'^, 
niB G {0,1}^^) into binary coded bits; ca = Ca (itia), 
Cb = Cb (hib), where C„ (•) is the binary channel code for 
user u, Cu is a length-L vector containing the coded bits of 
user u. The code rates of ca and Cb are va = EaIE and 
tb = LbIL, respectively. Next, we perform bit-wise XOR op¬ 
erations over the elements of ca and cb as the network coding 
operation; c^c = ca 0 cb. We then interleave the network- 
coded bits cnc using an interleaver tt; 'cnc = ’’"(catc), 
and lastly we map the interleaved bits Cjvc to the modulated 
signals. 

We assume the use of QAM, and its signal set is denoted by 
X. The size of X is \X\ = 2"*, where m is a positive integer. 
The signals in X are represented by complex numbers. The 
interleaved bits catc ure divided into N = L jm groups, with 
each group containing m bits. The i-th group is denoted by 

CNC,i = [cArc,(i-l)m-|-l: CAC,(i-l)m-l-2, ’ ’ ’ i ClVC.im] , where 
Cncj is the j-th element of catc, and i = 1,2, - ■ ■ ,N. Each 
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Fig. 7. The block diagram of our BICM scheme for rate-diverse wireless network coding. 


group of m interleaved bits is mapped to a modulated signal in 
X via a one-to-one labeling map p : {0,1}"* —> X. The i-th 
group Cjvc,i is mapped to the modulated signal x^c^i & X. 
We refer to ^ as the label of XMC,i- We use the same Gray 
labeling map /i as that used in single-user point-to-point BICM 

schemes | [T3] . The modulated signals XNC,i, * = !> 2, • • • , -/V, 

' ’ A T 

are stacked into a vector xjvc = [xnc,i,xnc, 2 , ■ ■ ■, xnc,n] 

for transmission to the users over the broadcast channel. 

The received signal at a users receiver is given by y„ = 
l^u^NC + Hu- We denote the j-th element of ync by yNC,i, 
the subset of the signals XMC,i G ^ where the j-th bit of 
the label of XNC,i (cNC,(i-^)m+j) has the value b G {0,1} 
by Xj^\ for j = and i = Based 

on yNc,i, the demodulator computes the detection metric of 


^ {yu,it C^C,(i—l)m+j — 

= Piyu,i\xNC,i)p{xNC,i\cNC,{i-l)m+j = b) 




(b) 


— Om-l 


I] P{yu,i\xNC,i) 


XNCyi^^. 


(b) 


(16) 

for b G {0,1}, j = 1, 2, • • •, TO, and i = 1, 2, • • •, TV. Then, 
these detection metrics are forwarded to the deinterleaver 
the output of which are the detection metrics of the network- 
coded bits {X{yu,i,C]s[C,i = (')}. where b G {0,1}, and I — 
1, 2, • • •, L is the Tth bit of catc- Using side information, users 
can easily deduce the detection metrics of the channel-coded 
bits of its target message. For example, user A obtains the 
detection metric of ca,i by using the side information cb,i as 


follows: 


(yA,i, ca,i = b) = X {yA,i^ cnc,i =b(B cb,i) (17) 

where b G {0,1}, I = 1,2, Finally, these detection 

metrics of channel-coded bits are input to the channel decoder 
for recovery of the target information bits. Some remarks about 
the above BICM scheme are in order: 

• The BICM scheme for rate-diverse wireless network 
coding adheres to the principle of virtual single-user 
channels expounded in Section III. Given the side infor¬ 
mation, the two users are actually served by two equiv¬ 
alent single-user point-to-point BICM schemes over two 
independent channels. Thus, single-user decoding after 
network decoding is sufficient to achieve the performance 
of the corresponding single-user point-to-point BICM 
scheme without losing optimality. 

• The data rates of users A and B are Ra = mrA, Rb = 
mrB- By employing different channel code rates ta, xb, 
the equivalent single-user point-to-point BICM schemes 
can flexibly support different data rates for user A and 
B. According to the point-to-point channel capacities of 
the users, Ca and Cb, we can tune the values of ta, xb 
and TO to optimize the system performance. In practice, 
the channel codes with different code rates ta, xb can 
be obtained from a mother code using methods such as 
puncturing |30| and splitting pT) . 

• Unlike our scheme here, |14| and (m made use of 
joint modulation for rate-diverse wireless network coding. 
Specifically, without considering channel coding, IE’ 
GD performed the network coding operations as joint 
modulations where the information bits of the two users 
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are jointly mapped to a constellation point. Using their 
different side information, the receivers of the two users 
can equivalently observe two sub-constellations with dif¬ 
ferent set sizes. Thus, these joint modulations use two dif¬ 
ferent constellations to deal with the different data rates of 
the two users. However, compared with the conventional 
constellations with the same cardinalities for single- 
user point-to-point channels, the minimum Euclidean 
distances between the points of the sub-constellations 
are shortened (See next section for details). This means 
that the equivalent single-user point-to-point SNRs are 
reduced. Therefore, the joint modulation schemes in 0. 
m do not conform to our principle of virtual single-user 
channels that guarantees two equivalent single-user point- 
to-point channels without any SNR reduction. By simply 
varying the code rates of the channel codes, our BICM 
scheme do not lose optimality, thanks to the principle of 
virtual single-user channels. The simulation results in the 
next section demonstrate the performance advantage of 
our BICM scheme over the joint modulation schemes in 
0 , 0 . 


VI. Numerical Results 

This section validates the performance of our schemes via 
simulations. 

A. The simulations of LDLC Rate-Diverse Wireless Network 
Coding 

This subsection presents simulation results on the perfor¬ 
mance of the LDLC implementation of the optimal frame¬ 
work for rate-diverse wireless network coding in section 
IV. In the simulation, the LDLC generating matrix with 
degree d = 7 is obtained from the generating sequence 

using the construction algo¬ 
rithm in 0. The codeword length is 100. Lor user A, the 
constellation of the elements of the integer message vectors is 
l^A,i = {~4, —3, • • •, 2, 3} for all i; thus the data rate of user 
A is Ra = log 2 I Li I =3 bits per channel use. Lor user B, the 
constellation of the elements of the integer message vectors is 
Lb,! = {~2, —1,0,1} for all j; thus the data rate of user B is 
Rb = log 2 I Li I =2 bits per channel use. The shaping method 
is the nested shaping method for LDLC. The BP decoding 
algorithm is used at the receivers of users A and B. Strictly 
speaking, the messages (the probability functions) in the BP 
decoding algorithm have real-valued variables. As in fTO) , we 
choose to use discrete vectors to approximate the messages. 
Specihcally, the i-th received signal is j/i = Xi -\- rii, the 
message associated with the variable node Xi is a continuous 
function fi{xi). We quantize Xi with a quantization step and 
over a range Xi G [j/i — A/2, ?/i-f A/2], where A is the 
magnitude of the range. In this manner, the message /i (xi) 
is represented by a vector. In our simulations, we use a 
quantization step 1/128 over a range of 8 in magnitude. Thus, 
each message is represented by a length-1024 vector. The BP 
decoding performs 100 iterations between the variable nodes 
and the check nodes of the LDLC factor graph to ensure 
convergence of the decoding process. 


We obtain the symbol error rates (SLR) of user A, 
P (bA,i ^ &A,i), and user B, p (bB,i 6 b.*), by averaging 
the decoding results over 1000 codewords for each user. The 
SER results of user A and user B in terms of SNR are shown 
in Lig. 1^ and Lig. respectively. Lrom the simulation results, 
we have the following observations: 

Lirst, the simulation results indicate that the LDLC imple¬ 
mentation of our framework for rate-diverse wireless network 
coding (RD-WNC) can indeed achieve the same performance 
as the LDLC for single-user point-to-point (SU-P2P) channels. 
Thus, as the codeword length goes to inhnite and the lattice 
generate by LDLC approaches the theoretically optimal lattice 
as discussed in Section III, the performances of LDLC SU- 
P2P will approach the Shannon limit fTO) , as well as that of 
LDLC RD-WNC. 

Second, consistent with our statement in Section IV, the 
shaping gain of LDLC RD-WNC (around 10 dB) is larger 
than that of LDLC SU-P2P (5 dB). 

Third, for user B in Lig. 11 (i.e., the user with the lower 
rate), LDLC RD-WNC has an additional lOlogm (6/5) « 0.8 
dB shaping gain compared with LDLC SU-P2P. This is also 
consistent with our analysis in Section IV.C. 

In conclusion, using implementable practical LDLC codes, 
we have empirically validated the feasibility and optimality 
of our framework for RD-WNC, which was shown to be 
theoretically optimal in Section III. 


B. The simulations of BICM Rate-Diverse Wireless Network 
Coding 

This subsection presents simulation results on the BLR 
performance of our BICM framework for rate-diverse wireless 
network coding (BICM RD-WCN) in Section V. In all simula¬ 
tions, the channel code used is the regular Repeat Accumulate 
(RA) code p^ . We can simply change the repetitions of the 
RA code to obtain different code rates. The decoding algorithm 
for RA codes is the iterative BP algorithm. The RA decoder 
performs 20 BP iterations to ensure convergence. The data 
packets for users A and B contain 10^ and 5 x 10^ binary 
source bits, respectively. Thus, the source rate for user A 
is twice that for user B. We compare our BICM framework 
with other schemes proposed for rate-diverse wireless network 
coding 0, 0. 

We hrst compare BICM RD-WCN with the joint modula¬ 
tion scheme in 0 - The joint modulation scheme originally 
proposed in 0 did not consider the effect of channel coding. 
In 0, the modulator of BS maps two bits targeted for user 
A and one bit targeted for user B (totally three bits) into 
a point of an 8-PSK constellation. The 8-PSK constellation 
labeling map is optimized to have large inter-point Euclidean 
distances at the receivers 0. The joint 8-PSK constellation 
and its labeling map proposed by |14| can be found in 
Lig.l (c) of 0. At the receiver of user A, using the 1- 
bit side information, the 8-PSK constellation is reduced to 
an equivalent QPSK constellation. At the receiver of user B, 
using the 2-bit side information, the 8-PSK constellation is 
reduced to an equivalent BPSK constellation. Although | [T4| 
focused on modulation and did not consider channel coding, it 
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SNR, 

Fig. 8. The SER results of LDLC RD-WCN and LDLC SU-P2P for user A. 



SNRb 

Fig. 9. The SER results of LDLC RD-WCN and LDLC SU-P2P for user B. 


is Straightforward to incorporate channel codes into the system. 
In our investigation here, for the joint modulation scheme, we 
apply an RA code with code rate r = 1/2 to the binary source 
bits of both user A and B. The coded bits are interleaved, and 
then the joint 8 -PSK modulation is performed according to the 
constellation labeling map shown in Fig.l (c) of GD- 

Unlike the joint modulation scheme, our BICM RD-WCN 
applies channel codes with different code rates to the two 
users, as illustrated in Fig. Specifically, RA codes with code 
rates ta = 1/2 and = 1/4 are used to encode the binary 
source bits of user A and user B, respectively. This results 
in the same total number of coded bits for users A and B. 
Bit-wise XOR are then applied to the coded bits of users A 
and B. The XORed bits are then interleaved, and the standard 
QPSK modulation with gray labeling map is applied on the 
interleaved XORed bits. Therefore, both BICM RD-WCN and 
the joint modulation scheme serve user A at the data rate of 
1 bit per channel use and user B at the data rate of half bit 
per channel use. The BER results of BICM RD-WCN and 
the joint modulation scheme are presented in Fig. 10 As the 
benchmarks, the performances of the SU-P2P BICM schemes 
(the r = 1/2,1/4 rate RA codes combined with the standard 
QPSK with gray labeling map) are also given in Fig. 

In Fig. 10 we can see that indeed BICM RD-WCN has 
the same performance as SU-P2P BICM; the joint modulation 


scheme and BICM RD-WCN have the same performance for 
user A; BICM RD-WCN gives better performance for user B 
than the joint modulation scheme does (around 3.5 dB gain at 
BER of 10-4). 

Eor the joint 8 -PSK modulation scheme, the equivalent 
QPSK at the receiver of user A is a standard QPSK constella¬ 
tion with gray labeling map. Therefore, it is easy to see that, 
for user A, the joint modulation scheme can achieve the same 
BER performance as BICM RD-WCN and SU-P2P BICM 
under the setup of the rate-1/2 RA code and the standard 
QPSK constellation. 

The 3.5 dB gain for user B obtained by BICM RD-WCN 
over the joint modulation scheme comes from two parts. The 
first part is the 0.7 dB gain from the modulation; the second 
part is the 2.8 dB gain from the channel coding. Eor the joint 
8 -PSK modulation, the constellation of the equivalent BPSK at 
the receiver of user B is not a standard BPSK constellation. It 
has an inter-point Euclidean distance d 2 = 1.85 (see in Eig. 1 
(c) of fM)). However, the inter-point Euclidean distance of 
the standard BPSK constellation is dspsK = 2. This implies 
that compared with the standard BPSK constellation, there is 
a lOlogj^Q « 0.7 dB performance loss in the 

equivalent BPSK constellation at the receiver of user B. It 
is well-known that the BER performances of standard BPSK 
and QPSK are the same in terms of Eh/No And, the 
modulation of our BICM framework is the standard QPSK. 
Thus, the first 0.7 dB gain in user B’s BER performance by 
BICM RD-WCN over the joint modulation scheme is from the 


modulation. Indeed, we can observe this 0.7 dB gain in Eig. 10 


(the gain of the SU-P2P BICM scheme with 1/2 rate RA code 
and QPSK modulation over the joint modulation scheme for 
user B). Our BICM RD-WCN accommodates the different data 
rates by adding more redundancies to the channel-cooed bits 
of the user with lower rate. In our setup, the channel code of 
user B is the rate-1/4 RA code for BICM RD-WCN and the 
rate-1/2 rate RA code for the joint modulation scheme. Prom 
the performances of SU-P2P BICM in Pig. 13, we know that 
the 1/4 rate RA code can provide about 2.8 dB more coding 
gain than 1/2 rate RA code does at BER of IQ-^. Therefore, 
compared with the joint modulation scheme, BICM RD-WCN 
obtains this second 2.8 dB gain from the channel coding. 

We next compare our BICM RD-WCN with the rate-diverse 
network coding (RDNC) scheme proposed in GD- Assuming 
the two users employ channel codes of the same rate, RDNC 
in fTS) just deals with rate diversity within the modulation 
process. We briefly review RDNC below. The binary source 
bits of user A and B are encoded into channel-coded bits 
using a channel code of rate r. The vectors of coded bits 
for user A and B are ca and cp. The length of ca is the 
double of that of cp. Then, the coded symbols are interleaved 
to obtain the interleaved bits: ca — 0 ( 0 ^ 4 ), Cg — n(cp). 
Then, to make the vectors of the interleaved bits have equal 
lengths, RDNC inserts zeros within cp in an alternate manner: 
cb,1j 0; cp_ 2 ) 0, cp_ 3 ,0, • • •, where Os are the inserted zero bits. 
The network coding of RDNC is the bit-wise XOR of ca 
and the zero-inserted Cp. Pinally, the network-coded bits are 
mapped to the constellation points of the modulation according 
to a special designed labeling map Gg. 
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Fig. 10. The BER results of our BICM RD-WCN and the joint modulation 
scheme IE) 



Fig. II. The BER results of our BICM RD-WCN and RDNC with QPSK 
in |15| . 



Fig. 12. The BER results of our BICM RD-WCN and RDNC with 16QAM 
in |15| . 


We investigate the BER performances of QPSK RDNC and 
QPSK BICM RD-WCN. The QPSK constellation labeling map 
for RDNC was developed in m (see Fig. 4 of (B)). Every 
two network-coded bits are mapped to one QPSK constellation 
point according to the specially designed constellation labeling 
map in i flS) . For QPSK RDNC, a rate-1/2 RA code is applied 
to the source messages of users A and B. For QPSK BICM 
RD-WCN, a rate-1/2 RA code is applied to user A, and a 


rate-1/4 RA code is applied to user B. The BER results are 
given in Fig. 11 We can see that, for both of users A and B, 
BICM RD-WCN gives better BER performance than RDNC 
does. For RDNC, the receiver of user A obtains a QPSK 
constellation after subtracting its side information; however, 
the QPSK constellation labeling map at the receiver is not 
gray labeling anymore. BICM RD-WCN does not lose the 
optimality of the SU-P2P BICM framework and the gray 
labeling for QPSK remains valid at the receivers. Due to this 
superiority in modulation, BICM RD-WCN gives better BER 
than QPSK RDNC does for user A, For RDNC, the receiver of 
user B obtains a standard BPSK constellation after substracting 
its side information; however, the channel coding of user B 
operates at a higher rate compared to BICM RD-WCN. Due 
to this superiority in channel coding, BICM RD-WCN also 
gives better BER than QPSK RDNC does for user B. 

We next consider the 16QAM modulation for BICM RD- 
WCN and RDNC. The 16QAM constellation labeling map for 
RDNC was developed in | [T5) (see Fig. 6 of |[T5)). Every four 
network-coded bits are mapped to one 16QAM constellation 
point according to the special designed constellation labeling 
map in d). The 16QAM constellation labeling for BICM 
RD-WCN is still the gray labeling. For 16QAM RDNC, a 
rate-1/2 RA code is applied to users A and B. For 16QAM 
BICM RD-WCN, a rate-1/2 RA code is applied to user A, and 
a rate-1/4 RA code is applied to user B. The above setups of 
16QAM RDNC and 16QAM BICM RD-WCN also give the 
same number of coded bits for both users. The BER results are 
given in Fig. 12 As the benchmarks, the performances of SU- 
P2P BICM schemes (the r = 1/2,1/4 RA codes with rates 
are combined with the gray labeling 16QAM) are also given 
in Fig. 12 We can see that our 16QAM BICM RD-WCN still 
has the same performance as 16QAM SU-P2P BICM; for both 
of user A and B, 16QAM BICM RD-WCN gives better BER 
performances than 16QAM RDNC does. The reasons are the 
same as the reasons for the comparison between QPSK BICM 
RD-WCN and QPSK RDNC. 

From the above simulation studies, we can conclude that our 
BICM RD-WCN can accommodate the different source rates 
of the users without losing the optimality of SU-P2P BICM for 
both users; other schemes proposed for rate-diverse scenarios 
cannot guarantee the optimality for both users simultaneously. 


VII. Conclusion 

We have proposed a nested-lattice-code encoding/decoding 
framework to achieve the optimal capacity pair of wireless 
broadcast channels with side information at the users. Al¬ 
though the nested-lattice-code framework is optimal theoret¬ 
ically, its exact implementation faces many difficulties. In 
particular, the lattice quantization decoding, which searches 
for the transmitted lattice point over the lattice space, has 
complexity that grows exponentially with codeword length 
N. Fortunately, with the principle of virtual single-user chan¬ 
nels suggested by our framework, we can implement our 
framework using practical codes with implementable decoding 
algorithms. To illustrate this, we implemented our encod¬ 
ing/decoding framework using FDFC that can be decoded 
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using a practical BP algorithm. We further applied the design 
principle to a BICM scheme, where binary channel codes 
are used. Simulation results of the LDLC and BICM im¬ 
plementations indicate that, for rate-diverse wireless network 
coding systems, our encoding/coding framework indeed can 
enable the two users to simultaneously achieve/approach their 
respective channel capacities. 
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